Stable topological textures in a classical 2D Heisenberg model 
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We show that stable localized topological soliton textures (skyrmions) with n2 topological charge 
f > 1 exist in a classical 2D Heisenberg model of a ferromagnet with uniaxial anisotropy. For this 
model the soliton exist only if the number of bound magnons exceeds some threshold value A^cr 
depending on f and the effective anisotropy constant KbS- We define soliton phase diagram as the 
dependence of threshold energies and bound magnons number on anisotropy constant. The phase 
boundary lines are monotonous for both v = 1 and v > 2, while the solitons with v = 2 reveal 
peculiar nonmonotonous behavior, determining the transition regime from low to high topological 
charges. In particular, the soliton energy per topological charge (topological energy density) achieves 
a minimum neither for 1^ = 1 nor high charges, but rather for intermediate values u = 2 or v = 3. 
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I. INTRODUCTION 

The studies of nonlinear excitations of two-dimensional 
(2D) and quasi 2D correlated spin systems are an im- 
portant issue of modern physics of magnetism, and 
can be useful for development of general soliton con- 
cepts lii^i^ii^i^ii The topological textures like localized soli- 
tons (skyrmions^) or magnetic vortices make an impor- 
tant contribution to the thermodynamics of magnetically 
ordered systems^ or even determine the character of its 
ordering as in the case of Berezinskii-Kosterlitz-Thouless 
transitioni^iiS Last years the interest for two-dimensional 
solitons has grown since they are frequently realized as 
ground state in the finite-size mesoscopic magnetic sam- 
ples, so-called magnetic dotsJ^ 

The most "famous" topological solitons are magnetic 
vortices having tti topological charge. These vortices 
are usually related to thermodynamic aspects of soliton 
physics i^ii^ Also, they appear in mesoscopic nanostruc- 
turesiii Easy-plane magnets with continuously degener- 
ated ground state have vortices with the energy being 
logarithmically divergent as a function of system size. 
The other and much less studied example of topologi- 
cal solitons are magnetic skyrmions which are present in 
isotropic or easy-axis magnets. Contrary to the above 
vortices, the latter textures are characterized by non- 
trivial 772 topological charge and finite energy. It is 
known that they determine the response functions of 
2D magnets at finite temperature a^^i^^ and take part in 
long-range order breakdown in isotropic mag nets^ The 
skyrmions form ground state of magnetic nanoparticles 
with easy-axis anisotropyii^ Their analysis is more com- 
plicated as compared to magnetic vortices and comprises 
many nontrivial features. An important example of latter 
features is the problem of a skyrmion stability since due 
to Hobart-Derrick theorem the static solitons with finite 
energy are unstable for wide class of models including 
standard continuous magnetic modelsi ^^'^^ 



For magnetic vortices, the consideration of lowest pos- 
sible topological charge = 1 is sufficiently as the vortex 



energy grows with 



Because of that it is 



advantageous for a vortex with i/ > 1 to decay for v vor- 
tices with 1^ = 1 and the vortices with = 2 can be stable 
in exceptional cases onlyJ^ The situation for skyrmions 
is not that simple. The simplest continuous model for 
isotropic 2D ferromagnet (FM) 
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admits the well-known Belvin-Polakov (BP) solutioufi^ 
which reads 
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where to is normalized magnetization 

TO = (sin 9 cos Lp; sin 9 sin Lp; cos ( 
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5 is a spin value, a is a lattice constant of a 2D FM. J is 
its exchange constant, r and x ^'I'e polar coordinates in 
the XF-plane, tpQ is an arbitrary constant. Solution ^ 
has the energy 

= i^Eo, Eo = 4ttJS^ (4) 

so that the state of BP skyrmions with v > 1 merges 
or dissociates into several other similar skyrmions with 
different v^s (the only rule that in such process the topo- 
logical charge should conserve) without their energy al- 
tering. Such exact degeneration is related to very high 
hidden symmetry, stemming from exact integrability of 
corresponding static model ([T]) (see, e.g.^^). This de- 
generation should certainly be lifted if we go beyond the 
model ([ij . The most important characteristic here is the 
parameter , which is appropriate to call the topological 
energy density, 

- — , (5) 
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where Ei, is the energy of a sohton with topological 
charge v. If S^, is a growing function of v, the most favor- 
able state with a given v comprises v solitons with unit 
topological charge otherwise such state is unstable. Lat- 
ter question is especially important for the investigation 
of general regularities of the highly excited magnet states 
evolution, (sec, e.g. Ref . [l9l and references therein) or for 
the analysis of essentially inhomogeneous magnet states 
under strong pumpingi^ Latter states can be generated 
by the ultrafast pulses of magnetic field, see Refs. [2lll2^ 
for details. The preceding discussion demonstrates that 
the problem of obtaining and investigation of the stable 
skyrmions with higher topological charges is extremely 
important.— 

The present work is devoted to the analysis of 
skyrmions with higher tt2 topological charges in 2D 
Heisenberg ferromagnet with uniaxial anisotropy ([6]). We 
show that there exists a certain range of system param- 
eters (exchange and anisotropy constants), where sta- 
ble precessional solitons with topological charge > 1 
exist. It turns out that in wide range of anisotropy 
constants, the topological energy density S^, of the tex- 
tures with ^ > 1 is lower then that of the textures with 
1/ = 1. On the other hand, the solitons with v = \ and 
V > 2 have monotonously growing phase boundary func- 
tions £i,,ci-{Ni,_ci-), while the case v = 2 has peculiar non- 
monotonous behavior, determining the transition regime 
from low to high topological charges. This means that 
the preferable values of soliton topological charge are nei- 
ther v = 1 nor high charges, but rather v = 2 or i/ = 3. 

II. MODEL DESCRIPTION AND SOLITON 
CLASSIFICATION 

We begin with the discrete model of a classical 2D 
FM with uniaxial anisotropy, described by the following 
Hamiltonian 

+A'^[(55)2 + (5|)2]. (6) 

n 

Here S = {S^ , , S^) is a classical spin vector with fixed 
length S on the site n oi a 2D square lattice. The sum- 
mations run over all lattice sites n and nearest-neighbors 
a, J > is the exchange integral and the constant k 
describes the anisotropy of spin interaction. In subse- 
quent discussion, wc refer to this type of anisotropy as 
exchange anisotropy (EA) . Additionally, we took into ac- 
count single-ion anisotropy (SIA) with constant K. We 
consider z-axis to be easy magnetization direction so that 
K > or K > 0. 

The analysis of real magnetic systems with discreet 
spins can be performed only numerically. In principle, 
it can be done by the same method as was described in 
Ref. [U, but it is not easy to extract necessary information 



from the set of numerical data. On the other hand, if we 
neglect the specific effects of discreteness (which appear 
at strong anisotropy only, K,k > J), like the presence of 
pure coUinear structures, see Ref. |24| . the consideration 
can be simplified using the generalized continuous ap- 
proximation and classical Landau-Lifshitz equation. In 
this approximation we can introduce the smooth func- 
tion S{x,y,t) instead of discreet variable Sn{t). In this 
case, the classical magnetic energy functional W[S] can 
be constructed expanding the discrete Hamiltonian ^ 
in power series of magnetization gradients, yielding 

^W2 + W4 + (7) 

where W2 contains zeroth and second order contributions 
to magnetic energy, see ([T]), and W4 contains the fourth 
powers of gradients. The explicit (and quite cumber- 
some) expressions for W2 and W4 had been presented in 
Ref. [2J. We note here, that single-ion anisotropy enters 
only W2, but not W4 and higher terms, while exchange 
anisotropy enters every term of the expansion ([7]) . 

For isotropic case K = and k = 0, W2 coincides with 
the energy of the isotropic continuous model ^ bear- 
ing BP soliton solutions of the form ^ with degenerate 
(with respect to topological charge) topological energy 
density (O, see also A simple accounting of mag- 
netic anisotropy in W2 generates a model, which is typi- 
cal example of the models governed by Hobbard-Derrick 
theorem - it does not admit static stable soliton solutions. 
In the model with W = W2 only, the size of any texture 
like domain wall, soliton etc (see, e.g. Ref. [Ufor details) 
is given by the characteristic length Iq 

ll = ^. K,n^K + 2t,. (8) 
2Kca 

In the case of weak anisotropy, K^s J, the length scale 
Zo 3> a so that the magnetization varies slowly in a space. 

Now we consider generalized model (O, including 
higher powers of gradients. In the expansion ([7]), we 
limit ourselves to the terms of fourth order only as they 
arc playing a decisive role in solitons stabilization, see 
Rcfs. j25i426l for details. The simplest possible general- 
ized model with account for W4 only can in principle 
admit the above static stable solitonsi^ Simple scaling 
arguments can demonstrate that. Namely, if the soli- 
ton core has size i?, the specific calculations with energy 
functionals like pB]) (see below) yield the dependence of 
the energy E^ on R in the form 

Eu = -A,^+B, + C,^, (9) 

where the first term comes from W4 and the rest come 
from W2- If E^{R) had a minimum (this occurs if 
Ai, < 0), this would mean the existence of a stable static 
skyrmion. Unfortunately, for the real model of magnet 
(ini) > so that the first term is negative and the 
dependence E^{R) does not have a minimum. We note 
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also that the expansion ([7]) and the above scale arguments 
are valid for any symmetry of initial 2D discrete lattice. 
The only difference is in the coefficients before gradient 
powers. Although these coefficients influence the solitons 
properties, the main feature of these expansions, consist- 
ing in the fact that in Eq. jO]) > remains the same. 
In other words, we did not find any symmetry of 2D FM 
with nearest-neighbor ferromagnetic interaction, where 
Ei,{R) has a minimum so that stable static soliton can 
exist 1^ 

In the absence of static two-dimensional solitons, it is 
possible to construct stable soliton states with station- 
ary dynamics, due to the presence of additional integrals 
of motion for magnetization fields. The purely uniaxial 
model ^ possesses the exact symmetry with respect to 
the spin rotation around z— axis so that the energy func- 
tional E,y[9, (fl does not depend explicitly on the variable 
ip. This leads to the appearance of an additional inte- 
gral of motion, z-projection of total spin. This integral 
of motion can be conveniently parameterized via integer 
N defining a number of bound magnons in a soliton N, 
see Ref . for details. In continuous approximation it can 
be written as 

TV = d^x{l-cose). (10) 

Conservation of N leads to the presence of so-called 
precessional solitons characterizing by time-independent 
projection of magnetization onto the easy z— axis and 
with the precession of the magnetization vector m at con- 
stant frequency uj around the z axis, 

= 9{r), ^ = u}t + iyx + ^o, (H) 

which holds instead of ([2|) in this case. The analogs of 
such precessional solitons are known to occur in differ- 
ent field-theoretical models, it is enough to note the non- 
topological Coleman's Q-balls^^ which do not have topo- 
logical properties as well as 7r2 topological Q - lumps )^ 
see Ref. d for details. 

Stable dynamical solutions with nonzero uj correspond 
to conditional (for fixed value) minimum of the energy 
functional . Namely, we may look for an extremum of 
the expression 

L = E^- fiujN, (12) 

where lo is an internal soliton precession frequency, which 
in this case can be regarded as Lagrange multiplier. Note, 
that functional is nothing but the Lagrangian of 2D 
FM magnetization field calculated with respect to spe- 
cific time dependence ([TT]) . This condition leads to the 
relation^ fiuj = dE^/dN , which determines the micro- 
scopic origin of the precessional frequency u. Namely, an 
addition of one extra spin deviation (bound magnon) to 
a soliton changes its energy by hu. Thus, the dependen- 
cies Ei,(N) and to^N) are very important for the problem 
of a soliton stability. 



III. THE STRUCTURE AND STABILITY OF 
SKYRMIONS. CRITICAL ENERGY. 

Further analysis of above continuous model consists 
in the solution of differential equations for soliton struc- 
ture, what is equivalent to the minimization of the func- 
tional ([7]). As these equations can barely be solved an- 
alytically, here we analyze the solitons properties by di- 
rect variational method. As we have shown earlier by 
comparison of variational approach and direct numerical 
minimization of initial discreet energy on a lattice)^ the 
variational approach gives fairly good results for weak 
anisotropics K^s < 0.5J, where the generalized continu- 
ous description is valid. 

The continuous models like ([1]) and ([7]) are usually 
parameterized by angular variables ([3]) so that the en- 
ergy E^ becomes a functional of these variables, Ei, = 
E,y[9,V6, (f>,W(j)]. Having energy functional i?^, we can 
write the corresponding Landau-Lifshitz equations and 
Lagrangian (fT^ . 

To apply direct variational method for minimization of 
the energy E^, = W2 + W4 we use the trial function 

Un- = ^j-^jJ-KAAr), (13) 

where K^{x) is the McDonald hmction with index ly. — 
Note, that trial function (fT^)) is based on the interpola- 
tive solution, constructed in Ref. |31i. The trial function 
(fT3|l gives correct asymptotics both for r (corre- 
sponding to BP soliton) and for r 00 (exponential de- 
cay with some characteristic scale 1/A), see Refs. 5,25,3l| 
for details. Latter exponential asymptotics is absent for 
solitons in isotropic 2D FM, e.g. for BP soliton 
It can be shown that the exponential asymptotics oc- 
curs for anisotropic models, where the length of decay is 
proportional to Iq ([8]). It had been demonstrated that 
due to power-law asymptotics of 9{r) for isotropic mag- 
net, the integrals defining soliton energy and number of 
bound magnons (see below Eqs. (jlSp ) are divergent for 
1/ = l .^-^^i'^^ To avoid this divergency in anisotropic mod- 
els like ^ , the interpolative (between BP asymptotics at 
r — > and exponential one at r — > 00) solution had been 
put forward in Ref. [3l|. Our analysis of continuous model 
([7]) shows that at v = 1 the above divergence can be cut 
off by the exponential asymptotics only. For 1/ = 2 the 
power-law asymptotics is sufficient for convergence of cor- 
responding integrals, while their derivatives with respect 
to parameter A are divergent. In this case, the diver- 
gences in derivatives are also eliminated by exponential 
asymptotics. At v > 2 all integrals and their derivatives 
(with respect to A) are convergent. This means that the 
behavior of solitons with v = I and 2 on one side and 
those with > 2 on the other side is different, being de- 
termined by the interplay between effects of anisotropy 
and higher spatial derivatives. 

In our minimization method, the parameter A is vari- 
ational, while the parameter R is kept constant as it is 
related to N, N (x R'^, see, e.gi^. In other words, we 
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minimize the energy E with the trial function (|13p over 
A for constant R. This approach has the advantage that 
it also permits to investigate the stability of obtained soli- 
ton texture. Namely, a soliton is stable if it corresponds 
to the conditional minimum of the energy at fixed A'^, and 
it is unstable otherwise. 

To proceed further, we introduce following dimension- 
less variables 

x = Ar, A = aA, z = AR. (14) 
and express trial function (fT5|) in terms of them. We have 

tan- = -^^j— yyy (15) 
I 



Then using the equation ((TO)) . we can calculate the num- 
ber of bound magnons in the soliton 

— = [1 - cos0)xdx = j^iIj{z). (16) 

In the variables p4|) the soliton energy assumes the form 

^ = ^7o(.)+7.W-^A^4(.), (17) 
where 



7o(z) = / sin 9xdx, 72(2:) = 7: xdx 
Jq 2 Jq 



5l'2 (l + Ksin^f 



zy^ sin^ I 



74(z) = / xdx\ {A^ef {l + K sin^ 6) +e'^{l + K cos^ 6) + 



sin^ ( 



x^ \x 



29'^ + 



x'^ 



dx ' dx'^ X dx 



(18) 



Thus, we express the energy and the number of magnons 
via two parameters, A and z. It turns out, that initial 
dimensional variables A and R enter the problem only 
in the form of their product z. The dependence of N 
and Ey on z enters the problem via a few complicated 
functions if), 70, 72 and 74, which can be written only im- 
plicitly in the form of integrals (|18p . However, in terms 
of these functions, the dependence on A (fT7|) turns out to 
be quite simple. This permits reformulation of the initial 
variational problem in terms of variables z and A^ only. 
Namely, we express 



A^ = 



2tt 



(N/S) 



4iz) 



(19) 



and substitute this expression into the dimensionless en- 
ergy ([T7)) . This gives us the expression for the energy 
of a soliton with given A'^, as a function of variational 
parameter z. Then we can find a minimum of E^ with 
respect to z, keeping A^ constant. 

The result of such numerical minimization in the form 
of the dependence £^{N) is shown in FigUfor a magnet 
with purely exchange anisotropy k ~ 0.05. The curves 
for higher k are qualitatively the same. To justify the 
applicability of our direct variational approach, earlier 
we have showu)^ that for 1/ = 1 the dependencies £1, (N) 
found by variational and numerical minimizations of the 
energy, are identical at small enough anisotropy, K^s ^ 
0.5J. 

Our analysis shows (see also FigllJ that topological 



(0 



K „=o.ij (k=o.o5J) v=i^ — : 




v=2 


v=3 



N/S 



FIG. 1: Dependence of a soliton topological energy den- 
sity £v{N) for skyrmions with different i/'s and exchange 
anisotropy only, k — 0.05J (T^eff = O.IJ). Dashed lines show 
£v,ci{i^) and N^^ci for — 1. 



energy density grows slowly as function of A^. We note 
here that this property holds for all anisotropy con- 
stants and topological charges. Even more interesting 
is the fact that all curves £u{N) have threshold points 
£u.cr ~ £u{Nu.ct:) SO that solitons exist only at £1, > £i,^cr 
and A"^ > 
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minimal soliton energy, which is most important char- 
acteristic of sohton contribution into magnet thermody- 
namics and can be observed experimentally.— For in- 
stance (see FigHI), at = 1 and k = 0.05 J f^^cr ~ 18 JS*^, 
which is substantially higher then that expected from Eq. 
Q 4ttJS^ « 12MJS^. 

Fig. [T]demonstrates one more unexpected soliton prop- 
erty, namely that at k = 0.05 J the energy density £v{N) 
decreases with increase of v and fixed N value. As we 
will show below, such behavior occurs for many values of 
i^eff, although there can be exceptions. The most impor- 
tant feature, however, is the existence of above threshold 
energy and bound magnons number values. The behav- 
ior of these values at variation of anisotropy constants 
is quite nontrivial. For example, at k = 0.05J (the 
value chosen for Fig. [T]) Ncr{v = 1) > Ncr{v = 2) 
but Ncr{v = 3) > Ncr{v = !)■ At the same time, 
the corresponding threshold energies behave monotoni- 
cally £ct:{^ = 1) > £ct:{^ = 2) > £cv{i^ = 3). Our ex- 
tensive analysis of numerical curves £v{N) for different 
anisotropics have shown that their overall behavior is dic- 
tated primarily by the threshold values: if Ey^ct > £v'.cn 
then the entire curve £y{N) lies above corresponding 
curve £u'{N) in wide interval of A^'s. This shows the 
importance of the above threshold values for the proper- 
ties of solitons. The information about these values can 
be conveniently represented in the form of so-called soli- 
tons phase diagram, i.e. the dependence of £„ and A^cr 
on AToff- 



IV. SOLITONS PHASE DIAGRAM. 

To obtain the phase diagram, we should pay attention 
to the details of above minimization procedure. Namely, 
to obtain truly stable soliton texture we should demand 
that the conditional extremum of the energy Ey is a min- 
imum. So, we keep track not only to the first derivative 
dEy/dz to be zero but also to the second derivative to 
be positive (corresponding to a minimum) at the point 
Zmin, where dE^/dz = 0. When a soliton approaches 
the limit of its stability, the modulus of second deriva- 
tive diminishes, becoming zero at the stability limit Zct, 
corresponding to above values £y = £y^cr and N = N^^cr- 
Thus the above soliton phase diagram is indeed deter- 
mined by the instability point (as a function of ratio 
Kcs/J) where both first and second derivatives of the 
energy are zero. 

The shape of the above phase diagram is determined 
by the character of anisotropy and v value and is reported 
on Fig.[2]for exchange anisotropy. First of all, one can see 
that the dependence £ct:{i') is quite complicated and its 
character change twice, at Kcs « 0.07 and at A'off ~ 0.29. 
It is seen from Figl5^ that while there is a quite large 
region of ATcff values, where £2 and £3 are smaller then 
£1 , the region where £3 < £2 is approximately three times 
smaller. Wc have shown that £4 is always larger then 
£3 although still smaller then £1. The energy £ct:{v) is 




0-^ ■ , ■ , ■ , ■ rJ 

0.0 0.1 0.2 0.3 I /0.4 



FIG. 2: Soliton phase diagram, (a) - threshold topologi- 
cal energy density (dashed line at small KcH is described by 
Ea. (|24|l ). (b) - threshold magnon number. Solitons exist at 
£ > £cr and A*' > A^'cr. For small k the value of A^cr is divergent 
as 1 / V Kes st V = 1 and 2 and is constant at v >3. 

growing with h' at K^s > 0.29 J only. 

Fig. [^b demonstrates the divergence of Ncr at small 
KcS as = 1 and 2. This is related to the fact that in BP 
soliton the integral describing A^ diverges logarithmically 
as r — > 00 for 1^ = 1, see also above. 

Overall, Fig. [2] demonstrates the complicated behav- 
ior of solitons phase diagram, which is needed to be 
understood. To understand better the above complex 
behavior, it is instructive to obtain the phase diagram 
analytically. Without loss of generality, we shall do so 
for the case of uniaxial anisotropy only. Such analytical 
treatment is possible since the anisotropy constant K in 
this case enters the problem only via coefficient in the 
first term of Eq. (|17p . This permits to obtain the ana- 
lytic dependencies £i,^cj:{K) and Ny^cr{K) hi the inverse 
form K{£,N). Note, that the structure of Eq. pB]) sug- 
gests that the above analytical procedure is applicable 
for both exchange and uniaxial anisotropy. Our analysis 
shows that this does not change the situation qualita- 
tively. Moreover, the quantitative results are close to 
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each other, compare Figs. [2]and|4l 

As it was shown above, the sohton phase diagram is 
determined by the instability point Zcr- To obtain an 
equation for latter point, we rewrite the Eq. (|17p in the 
form 



y = TT, a(z) = 7 ,/ 



(20) 



12 



^i{z)ij){z) 



(function is defined by Eq. (fT6|) ) and equate to zero 
the first and second derivatives of £y with respect to z. 
This gives the equation for the dependence y^Zcv) in the 
form 



2/c 



y(2cr) 



5'(z„)c"(Zer) 



C Uc 



.)6"(Zc 



a'(Zcr)c"(Zcr) " c' (Zcr )a" (Zcr) ' 



(21) 



where primes mean corresponding derivatives. Then, the 
equation for Zci can be obtained by substitution of above 
Uct: into the one of the equations determining zero of first 
or second energy derivatives with respect to z. It turns 
out that this equation always has real solutions if we 
formally admit the existence of negative K. 

The dependence (|2ip permits to obtain the equation 
for the phase diagram in the implicit form if (zcr), which 
reads 



K{z,,) 



J 



c"{z„) 



yl. + b"{z,,)y. 



2^{z) 



(22) 



The dependence Zci-{K), obtained by inversion of Eq. 
(im, is reported on Fig. [31 It is seen that while the 
entire curves Zcy{K) at z/ = 1 and 2 lie at K > 0, the 
curves for v > 2 \ic in this range only partially. Since for 
uniaxial anisotropy K can be only positive, only those 
parts of the curves with i> > 2 where K > correspond 
to physically realizable case. This clarifies the reason 
why the dependence Ncv{K/ J) for skyrmions with v > 2 
begins from finite N values, see Fig. [2b- This different 
behavior is related to the different kind of convergence 
of corresponding integrals dX v = 1 and 2 as well as for 
v > 2, see discussion above. 

The dependence Zci{K/J) can be easily recalculated 
to the solitons phase diagram £cii,{K/J) and Ncv{K/ J). 
This diagram is shown on Figs. [4^,b. It is seen the 
qualitative coincidence with the numerical curves from 
Figs[2^,b. The details of behavior can now be better 
seen then from above numerical curves. 

Our analysis shows that aforementioned different char- 
acter of convergence of integrals influences asymptotics 
of the phase diagram curves at small anisotropics. This 
influence can be seen from Figs. [5^ and [H where the 
dependence S^^criK/J) at v = I and 2 is nonanalytical 




FIG. 3: Dependence Zci{K/ J). Numbers near the curves cor- 
respond to v values. 



at [K/J] — > , while for > 2 it is analytical. The 
reason for such behavior can be seen from Fig. [Sj show- 
ing that small K/J correspond to small z for v = 1 and 
2, while ior V > 2 all values of K/J including limit- 
ing case K/J correspond to flnite z. This means 
that at I' > 2 the asymptotic analysis of the curves 
£v,ci{K/J) can be done simply by Taylor expansion at 
small K/J, which yields simple monotonic behavior of 

^f.cr and Nu.cr- 

At = 1 and 2 the situation is not so simple and re- 
quires more complicated analysis. Such analysis can be 
performed on the base of Eqs. (HI]), (|^ and is quite 
cumbersome. The main idea is that functions a(z), b{z) 
and c{z) in (|20p at small z can be represented via ex- 
ponential integral functionsi^ which, in turn, may be 
expanded in asymptotic series. This procedure for = 1 
gives following parametric dependence 



K _ 18 

y ^ hTs 

iVcr _ 47r 
S 



Tn^z, 



1 



Sc. 



AttJS^ 



1 



In 3 



'In^ 



9 z^ln' 



(23) 



The asymptotics is shown by dashed line on 

Fig.[4li. In the main logarithmic approximation this gives 
simple relations 



£i,cr = ^ttJS^ [l + Ci VK/Jj , (24) 

SC2 



(25) 



which coincide with the results obtained in Ref. |24| by 
direct numerical modeling of the discreet model, giving 
Ci « 1.87 and C2 ~ 5.65. 

These expressions give us an idea about behavior of 
energy and bound magnons number at small K. Namely 
we see that the energy has the (approximate) square 
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cases V = \ and = 2, which can also be seen from 
Fig. m Namely, if the topological energy densities for 
skyrmions with v ^ \ and 2 have both square-root non- 
analyticity, for v = 1 this nonanalyticity reveals only as a 
small cusp in a close vicinity of if = 0. Other peculiarity 
of = 2 case is nonmonotonous behavior of £c-ci (decreas- 
ing at K < O.IJ with subsequent increase). Our analysis 
shows, that the above peculiarities of the skyrmion with 
1^ — 2 arc due to interplay between exponential asymp- 
totics of 0{r), necessary to cut-off the divergence at v = 1 
and 2 as well as power-law asymptotics, which is sufficient 
for the convergence a.t v > 2. 
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FIG. 4: Soliton phase diagram for the model with uniaxial 
anisotropy only, (a) - threshold topological energy density 
(dashed line at small Kj J is asymptotics p3|l 'l. Inset expands 
the encircled area on the main panel and shows the behavior 
of £2,cr at small Kj J . (b) - threshold magnons number. The 
ranges where solitons exist are similar to those on Fig. [2] 
Numbers near curves correspond to v values. 

root singularity, while magnons number A'^ obeys in- 
verse square root law. The asymptotics ([^^ is shown 
on Fig. [4^ by dashed line. 

For V ~2 the same approximation gives 

K _ 0.0066482 ^ 0.0132207z 

J \y? z \\? z 

^ = -187.047 Inz - 108.007, 

£cr , , 0.02762 , 0.01595 
4^ = ^ + n^^ + 17^- 

This also yields the divergent N2^cr{K) and square root 
peculiarity in the energy £^2,cr(-^)- However, the coeffi- 
cient before y/K/J in the energy is much smaller then 
that a,t V = 1 , which makes it almost invisible in the scale 
of Fig.dl The details of this behavior are reported on the 
inset to Fig. 2^. 

This shows the similarities and differences between 



V. DISCUSSION AND CONCLUDING 
REMARKS 

In this paper we present a comprehensive theoretical 
study of the localized topological solitons (skyrmions), 
stabilized by precessional spin dynamics, for the clas- 
sical 2D ferromagnet with easy-axial anisotropy on a 
square lattice. Our efforts were directed primarily to- 
wards the study of the role of higher 7T2 topological 
charges v > 1 on the above solitons properties. Our 
main conclusion is that the interplay between high topo- 
logical charges, effects of lattice discreteness (in the form 
of higher powers of magnetization gradients) and uniaxial 
magnetic anisotropy makes many unexpected peculiari- 
ties into solitons properties as compared to those in the 
simplest isotropic continuous model with W4 = 0, con- 
taining BP solitons. 

Similar to previous studies, it turns out that the pres- 
ence of even weak anisotropy makes solitons dynamic, 
i.e., those with nonzero precession frequency for any num- 
ber N of bound magnons. The minimal consideration 
of discreteness (via higher degrees of gradients of mag- 
netization) yields the existence of some threshold value 
of both soliton energy (topological energy density (O, 
which is more appropriate characteristic for solitons with 
h' > 1) and the number of bound magnons. Similar to 
the problem of cone state vortices^^ the instability is re- 
lated to the joint action of discreteness and anisotropy. 
As a result, the critical values of bound magnons TV^.cr 
and soliton energy £i,^cr is present and nonanalytic de- 
pendencies of above threshold values on the anisotropy 
constant appear at ly ~ I and 1/ = 2. It was shown 
earlierj^l that the variational minimization of W2 (i.e. 
energy, incorporating only squares of magnetization gra- 
dients) gives no threshold for soliton existence, i.e. the 
soliton exists everywhere up to TV = 0. This means, 
that mapping of the initial discrete model even for small 
anisotropy A'ofr ^ J on the simplest continuum model 
(HI) is wrong and to get the correct description of solitons 
in 2D FM we have to consider at least fourth powers of 
magnetization gradients. This seemingly paradoxical re- 
sult is actually due to the fact that the terms with (Vm)^ 
are scale-invariant (so that the corresponding energy has 
a saddle point) while the (stable) soliton size is deter- 
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mined by the fourth derivatives as well as by a magnetic 
anisotropy. Our analysis show, that higher (then fourth) 
powers of magnetization gradients do not change the sit- 
uation qualitatively, rather, in the range of above-studied 
K^ff these terms make only a small quantitative contri- 
bution to the solitons phase diagram. This means, that 
solitons with = 1 can be well studied within the model 

0- 

In summary, we have shown the existence of stable 
topological TT2 solitons (skyrmions) in the 2D ferromag- 
net with uniaxial anisotropy on a square lattice. Since 
for any 2D lattice symmetry the structure of correspond- 
ing energy functional ([9]) is the same, we may speculate 
that such textures exist in any 2D magnet with uniaxial 
anisotropy. We also show that due to different character 



of the energy and magnons number integrals convergence 
at 1/ < 2 and i/ > 2 determines the features of so-called 
solitons phase diagram. The main nontrivial and un- 
expected result of the paper is that while the solitons 
with v = 1 and ly > 2 have monotonously growing phase 
boundary functions £^i/,cr(-/V,y_cr), the case ly = 2 has pecu- 
liar nonmonotonous behavior, determining the transition 
regime from low to high topological charges. This means 
that the designated value of soliton topological charge, 
which is expected for highly excited state of FM, is nei- 
ther V = I nor high charges, but rather = 2 or i/ = 3. 
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